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A mechanism for the formation of a new type of stationary state with macroscopical number of 
phonons in condensed atomic gases is proposed. This mechanism is based on generating longitudinal 
phonons as a resuh of parametric resonance caused by a permanent modulation of the transverse trap 
frequency in an elongated trap. The phonon-phonon interaction predetermines the self-consistent 
evolution which is completed with macroscopic population of one from all levels within the energy 
interval of parametric amplification. This level proves to be shifted to the edge of this interval. All 
other levels end the evolution with zero population. 
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I. INTRODUCTION 

At zero temperature the formation of the stationary 
state in an atomic condensate with involving of macro- 
scopical number of phonons certainly requires thermody- 
namically non-equilibrium conditions. 

In the present paper an idea is proposed for realizing 
these conditions for longitudinal sound phonons in elon- 
gated Bose-condensed ultracold gases. The idea is based 
on generating phonons as a result of the parametric reso- 
nance induced by a permanent modulation of the trans- 
verse trap frequency uj_\_. In last years, it has become 
clear that the parametric resonance plays an important 
role in a number of phenomena in ultracold gases. As is 
shown in [l[, the parametric resonance, responsible for 
the effective energy transfer between different branches 
of excitations, determines the temporal evolution of the 
system, in particular, origin of the damping for the trans- 
verse breathing mode at zero temperature. The further 
consideration of the nonlinear evolution accompanying 
the parametric resonance [S] allows to explain the non- 
trivial picture of the damping revealed in the remarkable 
work of the Paris group [3[ . The parametric amplification 
of the Bogoliubov modes due to periodic modulation of 
an optical ID lattice is considered in [1]. The interesting 
picture of the density distribution in low-dimensional [5| 
and toroidal [ij] Bose-condensed gases under conditions 
of the parametric excitation has been revealed at the nu- 
merical solution of the Gross-Pitaevskii equation. 

A permanent modulation of the transverse trap fre- 
quency Ci;_L brings about the oscillations of the conden- 
sate density and, accordingly, sound velocity. Owing to 
the parametric resonance, this results in the generation 
of pairs of sound phonons with the opposite momenta 
and energies close to a half of the modulation frequency 
loq /2. In fact, phonon pairs are produced within the finite 
energy interval a;o/2 ± Sq. In what follows, this inter- 
val of parametric amplification with the width 2Eq near 
wq/S is identified as "the parametric interval". The pa- 
rameter Eq is connected with the modulation amplitude 



|(5aj^/a;^| <C 1 by the relation Eq ^ ujq ■ \Suj±/uj±\. If Eq 
exceeds the damping factor 7 for longitudinal phonons, 
their occupation numbers begin to increase exponentially. 
In this case the phonon-phonon interaction becomes sig- 
nificant. As a result, the evolution scenario changes dras- 
tically. 

For the typical longitudinal size of traps L a few 
phonon levels fall into the narrow energy interval 2Eo. 
The main effect of the phonon-phonon interaction, in 
addition to the damping, manifests itself in an effec- 
tive renormalization of the levels. In the course of the 
self-consistent nonlinear evolution the position of lev- 
els changes continuously with changing the number of 
phonons. As is shown below, in the simplest case of one 
doubly-degenerate level in the parametric interval 2Eq, 
the evolution is over as the renormalized level reaches the 
left or right boundary of the interval. Then the growth of 
the phonon occupation numbers ceases. As a result, the 
steady state of phonons with the macroscopic occupation 
number in a single state and the common phase for all 
excitations appears. 

In the case of several levels the picture is more sophis- 
ticated. It turns out that only one of the phonon levels 
finishes the evolution achieving the stationary state with 
the macroscopic number of phonons. As in the case of 
a single level, this level proves to be at the edge of the 
parametric interval. 

The preliminary analysis of the problem for the case 
of a single level within the parametric interval has been 
represent in our work [8]. In the present paper we give 
a comprehensive consideration of the problem, including 
the general analysis of the evolution for arbitrary number 
of phonon levels within the parametric interval. There 
is considered the essential problem of indrawning outer 
levels located initially out of the energy interval of para- 
metric amplification. Since in general case the evolution 
is over with the stationary occupation only a single level, 
the question about stability of this state is specially anal- 
ysed. 
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II. PARAMETRIC GENERATION OF PHONON 
PAIRS 

1. Let us consider the Bose-condensed gas at T = in 
a trap of the cyHndrical symmetry with L 3> i?, where 
R is the radius of the condensate. Neglecting the edge 
effects, we can write the general equation for the field 
operator ^'(r, z, t) of atoms in the form 



ih 



5* 

dt 



2m 



2m dz^ 



■^+Uo-^+M, (1) 



where Uq — Anfi^a/m, a being the 3D s-scattering length. 
Note, considering a dilute and cold Bose gas with domi- 
nant binary collisions wc imply as usually that the radius 
of interaction region 7'o ^ a, the gas parameter na^ <C 1, 
and the correlation length ^ > n~^^^ {n is the atomic 
density). In this case, the replacement of two-body in- 
teraction with the effective contact potential is relevant. 
We assume that the trap transverse frequency depends 
on time as 



UJ±{t) — UJ ±{1 + rj sin ujQt), 77 <C 1. 



(2) 



According to Ref. |8|, we introduce the spatial scaling 
parameter b{t) and new variables ~p = r/b{t), T{t). The 
field operator in terms of new variables can be written as 



bit) " ' 



2h 6' 



(3) 



Substituting these expressions into Eq.(IT]), we obtain 
straightforwardly 



OT 



2m^P 



2 , muj^p^ 



2m dz^ ' 
(4) 

provided that functions b{t) and T{t) satisfy the equations 



where /z is the chemical potential corresponding to the 
static cylindrical trap with wj, = const. In the Thomas- 
Fermi limit when nU^ ^ fujj±^ the solution of this equa- 
tion is reduced to known expression 



xl 



Uo 



1 



2fi 



In the quasi-lD case when nUo <^ hLU±, we have 



xl 



ni 



■ cxp 



ilh^{t)r 



mijj±^ 



In the both cases the modulation of the trap transverse 
frequency results in the oscillations of the condensate 
density. 

2. Considering excited states, we first linearize Eq. (j4|) 
in x'. Then from Eq. (4) we have 



dr 



2m dz^ 



ho 



V 

2m P 



muj\(P' 



(8) 



Under condition <C w_l only long-wave longitudinal 
phonons prove to be involved into the evolution of the 
system. 

For these phonons the condition kR <C 1 is fulfilled and 
the transverse distribution is close to xo (sec, e.g., [9]). 
Comparison of Eqs. (7) and (8) allows us to simplify 
the latter equation omitting the first term in the right 
side. Using the second relation in ^ and going over to 
variable we arrive at the following equation 



dt ~ b-'ity^ '^^ ' 2m dz^ 



(9) 



Involving the symmetry of the problem, we have for the 
operator x' in second quantization 



fb 



ujl{t)b 



Ul 



± . 

53 ' 



, dr 
'dt 



1. 



(5) 



Treating the condensate wave function as independent of 
z, we find that Eq.Q being expressed in variables p, r de- 
scribes the evolution at the time-independent frequency 
Lu±. Combining the expression ([2|) with the condition 
Wq <C w± , we find the solution of Eq. ([5]) in the form 



b{t)^l + bi(t), bi{t) 



sincjo^- 



(6) 



The field operator x can be represented in the usual form 
X = (Xo + x') exp(— z/Ltr), xo being the condensate wave 
function and fj, is the initial chemical potential. In the 
absence of excitations at T = we can omit the last term 
in Eq. ([4]) and determine Xoip) from the equation 



, Xfeflfe; Xfc 



^•</>o(7?), (10) 



where being the annihilation operator for atoms. At 
ujQ <C Lu± the function (j)o is actually close to xo for both 
the quasi- ID and Thomas- Fermi cases. In the expression 
PU)) we keep the ground state alone for the transverse 
motion assuming that higher states are insignificant in 
the course of evolution. Substitution of the expression 
into Eq.® gives after the standard averaging over 
the radial variables 



ih- 



dcLi^ 

~dr 



G- 



2m 



ak+Gatk-vG sin LJot{dk+at^.) (11) 



— V - 
2m P 



G = f d^ p(I)q{~p)G{~p). Let us rewrite this equation in 
terms of phonon operators using the usual transforma- 
tion dk — Ukbk + Vko^k- When the last term in Eq. (jlip 
Xo + C^Xo — 0' (''') is omitted, a set of equations for 6^,6^1. determines the 
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well-known Bogoliubov spectrum ojk (with the replace- 
ment nUo G). Considering the nonstationary case, we 

introduce the substitution hk = hk exp(—iuJot/2) and take 
into account that long times loqI 3> 1 are most interesting 
for the analysis. Then, within the resonance approxima- 
tion, which implies that we discard the strongly oscillat- 
ing terms with exp(±iLJoV2): the equations for 6fc,fe^^ 
are reduced to the form 

dbk ~ -I "^^ifc ~+ - / N 

i-j^ ^ S,kbk+iEokbZk\ -i—^ ^ £.kbZ,^-iEakbk, (12) 

Eak = {i-jG /2h)[{uk + Vk) /uk], Cfe = For sound 

phonons |(ufc + Wfe)/wfe| « huk/n- Since Eok ^ ojq and we 
are interested in the energy interval jwfc ^ luo/2\ < 2Eo, 
the frequency uJk in the expression for Eok can be replaced 
by ujQ. In the other words, Eok is replaced by its value 
at k = |fco|, where |fco| is determined from the condition 
ujQ = 2u}ko- ^ result, we have 

\6uj I 

Eoko ~ 11^0 ^ i^o = Eq. (13) 

Substituting the transformation bk = 
cifc exp(zeofct)6fe + C2ke:x.p{~ieokt)btk into (HH), we 
find EQk = iaok with aofe = [|£^oP - S.IV^'^- This 
result demonstrates the appearance of the parametric 
resonance with an exponential growth of the phonon 
occupation numbers at \Eo\ > which is induced 

by the modulation of the transverse trap frequency. 
Thus the parametric resonance occurs within the narrow 
range near luq with the width 2Eo. 

III. EVOLUTION EQUATIONS 

1. So far we have disregarded the phonon- phonon in- 
teraction. With an exponential growth of the phonon 
number the interaction begins to play an essential role. 
The weakness of the phonon-phonon interaction as itself 
(see below) allows us to take only three- and four-phonon 
processes into account. Assuming generation of sound 
phonons alone in the parametric interval 2Eq, we can use 
the expressions for iJ*^^^ and iJ'*'_obtained within the hy- 
drodynamic approximation, see [ll|. Owing to smallness 
of the gas parameter, the dominant term in the Hamil- 
tonian of three-phonon interactions has the form [ll| 

Here v, 5n are the operators of velocity and alternating 
part of the density, respectively. In fact, we consider the 
interaction of longitudinal sound phonons that effectively 
reduces H'^^^ to the one-dimensional problem. 

Using the known expressions for operators v, Sfi [ll[, 
one can find the vertex A3 for i?^^^ in (14). 



Since each three-phonon vertex implies the momen- 
tum conservation law, at least one of three phonons lies 
beyond the parametric interval and, therefore, has zero 
occupation number at T = 0. This makes possible to re- 
duce the expression obtained in second order in H^^^ to 
the effective Hamiltonian for the four-phonon interaction 
(A(fc) is the Kroneker symbol) 

= 4 ^ b+^blbk,bkAiki + k2-h-h). (15) 

Here all states lie within the interval of 2Eo. Using the 
expression for A3 and taking into account that all wave 
vectors \ki\ (Eq ^ uq), we find (to the accuracy of 

numerical factor of the order of unity) 

|A| . ^. («) 

where N is the total number of particles. 

The same four-phonon interaction can be obtained as 
a result of the standard canonical transformation of the 
Hamiltonian Hq + H'^^'> (see, for example, [l^). 

The direct calculation shows that the ratio of the ma- 
trix element in H^"^^ from Ref. [ll| to \A\ is proportional 
to (na3)i/2 < 1. 

In addition, second order in 7?'^' contains the imag- 
inary part related to real decay processes of phonons, 
which determine the phonon damping. The later we take 
into account phenomenologically introducing a decre- 
ment 7. After work [3| one can conclude that the pa- 
rameter 7 is small at T = for the geometry under 
consideration. Therefore, it is rather easy to satisfy the 
conditions when aofe > 7 and the parametric growth of 
the number of phonons remains. 

2. As is follows from the above expressions, the Hamil- 
tonian of interacting phonons near the parametric reso- 
nance have the form 

H = J2^kbtbk+J2{Eoe-'-°'btbtk + H.c.)+H,ff (17) 

k k 

Let us find (dbk/dt) taking into account the phonon- 
phonon interaction. Using the substitution bk ~ 

bk exp(— iajoi/2)i we have 

=(-*7/c + ^k)bk + iEobtk+ 

^ ^ ^ (18) 
-I- Yj AblbkJ)kA{k + k2~'k3-ki). 

k2k3k4 

By means of the equations for bk and b^f^ one can directly 

obtain the equations for the correlators Nk b'j^bk > 

and fk —< bkb-k >■ These equations contain the four- 
phonon correlators. Within the mean-field approxima- 
tion we represent the four-phonon terms as products of 
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the two-phonon correlators Nk and fk- This imphes that 
we keep the terms linear in A. In addition, we go over 
to the classical Bose field for phonons supposing that 
the number of parametrically excited phonons is large. 
As a result, we arrive at a set of nonlinear equations 
that describes the self-consistent evolution of interacting 
phonons within the parametric interval. These equations 
have the form 



dNk 
dt 



^=-2{j+t^k)fk + 2EoNk'2iAVNk, 



(19) 



where the position of the renormalized level is determined 
by the expression 



^ + Q + ANk; A^3A. 



(20) 



Here is the level position at the initial time moment, 
Q^2A J2 Nk';r= J2 /fe'- Hereafter the evident 

k'^±k k'^±k 

relations ujk = ^-k, Nk = N^k are used. 

For the large occupation numbers Nk to the accuracy 
of terms ~ ^/Nk, the function fk{t) can be written in the 
form fk{t) = \fk{t)\exp[i(pk{t)] = Nk{t) exp[i(pk{t)]. As 
a result, we re-arrange (fTO|) to a set of equations for Nk, 

Vk 



dNk 
dt 



= -2-fNk+2EoNkCOSipk- 



2ANk Nk' sm{^k - fk'); 

(21) 



d(pk 
dt 



= —2^k — 2E(jSinLpk-~ Nk' cos{ipk 

k'^±k 



fk' 



Hereafter we take into consideration that ipk = ip-k- 

In what follows, we consider the case of finite longi- 
tudinal size L assuming that a single or a few discrete 
two-fold degenerate levels lie within the energy interval 
of about 2Eo. 



IV. SOLUTION FOR A SINGLE LEVEL 

1. In the case of a single level the equations (pij) take 
the form 



dNk 
dt 

d(pk 
dt 



-2jNk + 2EoNkCosipk; 



(22) 



-2[^k + ANk + Eosmipk]. 
This system has the following stationary solution 



VW^T^k 



; sm ifi^ = T 



En 



(23) 



The signs =F correspond to ^4 ^ 0, respectively. For def- 
initeness, we suppose Eq > 0. The result has an 



interesting physical origin. In fact, the interaction -ffe// 
determines the effective renormalization of the phonon 
level, leading to Suik — ANk for the case concerned (see 
Accordingly, ^ = ^k + ANk- From Eqs. 
(Pg)) it straightforwardly follows that = ±^/E^~^. 
So, at Nk — N^ the renormalized level reaches the left or 
right edge of the parametric energy interval (within the 
accuracy of the shift due to 7). As a result, the para- 
metric increase of the phonon occupation numbers stops 
and the phase acquires the constant value. The maximal 
value of the phonon number equals 



Nt 



Eo 



7,^fe < Eo 



(24) 



The phase is 93^ ~ T'^/2 for A ^ 0. The phase ipk corre- 
sponds to the phase correlation of phonon pairs with the 
opposite momenta. The appearance of the phase ipk and 
correlator fk is connected, first of all, with the creation 
of phonon pairs in the course of the evolution as a re- 
sult of the parametric resonance. The evolution finishes 
in the stationary state of phonon pairs with zero total 
momentum and the common phase ip'l.. 

2. The character of the temporal evolution, tending 
asymptotically to the values (^5]) . ([M)) . depends essen- 
tially on the relations between the parameters in Eqs. 
([22)) . First, let 7 = 0. In this case, Eqs. ^ have the 
conserved integral of motion 



Ho = 2EoNk sin (pk 



ANl 



2^kNk 



dHo 
dt 



0. (25) 



The variables Nk,fk are canonically conjugate, and the 
equation of motion (|22p can be written in the Hamilto- 
nian form 



dNk 
dt 



SHq 

Sfk ' 



dipk 
dt 



SHq 

SNk 



At the initial time moment, when the modulation of uj±_ is 
switched on, we suppose ^fc(O) = £,k + ANk{0) and aok = 
\jEl - ffc(0)2. At that, Eqs.dUl) imply that cos¥)fc(0) = 
(iok/Eo- Bearing in mind that sin(^fc(0) — — ffe(0)/£^Oi 
from the second equation in ([22]) we obtain {d(pk / dt) (0) = 
0. As a consequence, the initial (and conserved) value of 
Ho is equal to 



Ho ~ -AN^iO) < Eo. 



(26) 



For an arbitrary time, Eq. (|25p can be rewritten as 
2Eosmipk{t) = -[2Ck + ANkit)] - Ho/Nk{t). Substi- 
tuting this relation into equations (|22p . one can find 



dNk 
dt 



±2Nka; 



d(pk 
dt 



-ANk 



Ho 

Nk' 



(27) 



a = [E^ - {ANk/2 + a - i/o/2iVfe)2]i/2. The signs ± 
correspond to the regions with \ipk\ < 7r/2 and l^s^l > 



7r/2, respectively. For the region with \Lpk\ < 7r/2, the 
solution can be found straightforwardly 



2t: 



dx 



Nk{0) 



[i?2-(Ax/2 + 6)2]V2' 



(28) 



Here we omitted the small term with \Ho\ ^ Eq. The 
upper limit of the integral is equal to the value N™ w 
2{Eo =F at which the denominator vanishes. The 

divergence at Nk ^ is a typical manifestation of the 
parametric resonance which requires a nonzero initial 
field amplitude. This can be achieved by taking zero- 
point oscillations into account fi\. Owing to these oscil- 
lations, we can put Nk{0) ^ 1. The time necessary for 
the system to achieve the maximal value N™ is equal to 



1 



2a 



In 



Ofc 



A 



^2 
^0 



(29) 



We see that the argument of logarithm is much greater 
than unity at < Eq and Eq ^ A. This implies that 
tm S> l/2aofc where l/2aoA: is the characteristic time of 
the parametric resonance. At t = tm we have \ipk\ = 7r/2. 
As it follows from the second equation in ((27l) . a,t t > tm 
the phase proves to be in the region with \ipk\ > n/2. 
As a result, in the first equation the sign becomes neg- 
ative and the phonon number reduces. This is the start 
of an oscillating behaviour. The numerical simulation 
shows that the joint solution of Eqs. ((27)) has the form of 
anharmonic oscillations of both the phonon number and 
the phase around their stationary values. At £^k = the 
exact solution for Nk{t) can be written in terms of dn 
Jacobian elliptic function as 



Nk (t) = NJ!'dn[2Eo {t-tm);k], NJ^ 



2Eo 



(30) 



Where dn{u\ k) = dn(—u] k), dn{0; k) — 1; k is the ellip- 
tic modulus (l3| . For the given initial conditions we have 

yr^p = \hqA\/ae^ < 1. 

3. The character of evolution changes drastically at 
7 7^ 0. Now in all cases the system asymptotically ap- 
proaches the stationary state. At any set of initial param- 
eters the solution of Eq. (|22p leads to the values (|24p for 
these stationary states. The damping of the oscillations 
has the decrement, which is close to 7, when ao S> 7 and 
t ^ tm • A.S 7 increases, the arrival time to the stationary 
state decreases as I/7. When 7 is relatively close to aok, 
the oscillations disappear. In this case the time necessary 
to reach the stationary state is allied to the value (p9)l . if 
aok is replaced by ai = aok — 7. 

The direct numerical simulation of the system ((22|) 
demonstrates the picture described. As an example, in 
Fig.l we report the results for the phonon number Nk 
as a function of the dimensionless parameter s = 2jt 




FIG. 1: Phonon numbers Nk{s) versus s — 2'yt (a,b) and 
versus y{t) = ~^p{t) (c,d) for (£0/7) = 5, 1.3. 



for E^l-i = 5 (a) and EoJ-f = 1.3 (b). Figs. 2 c), 
d) demostrate the phase portraits in the plane [Nk^i^k] 
for these cases. (The parameters {Eq/A) = 10^ and 
^ = are used for the calculations.) In the both cases 
Nk{t), ipk(t) are assymptotically approaching to the val- 
ues 

Thus at any finite values of 7 the self-consistent evo- 
lution of interacting phonons experienced the parametric 
amplification results in the formation of the stationary 
state with the macroscopic number (of the order of Eq / A) 
of phonons in a single state. 

It should be noted that the character of the temporal 
evolution is sign-independent of A. 

It is interesting to verify the stability of the solution 
obtained. It seems that the damping of fluctuations of 
the phonon number SN = Nk — N^ may have a differ- 
ent character within and beyond the parametric interval. 
Really, at SN < the phonon number reverts to the sta- 
tionary state owing to parametric pumping. While at 
SN > 0, the absorption alone, determined by 7, returns 
the phonon system to this state. 

Linearizing Eqs. ()22p in SN and Sip ~ ipk — iPk with the 
initial conditions SN{0) = Nk{0)~N^ and Sip{0) =^ 0, we 
obtain 

^'^^ -[2EoN^sm^l]Sip; 

(31) 



dt 



^ = ~2ASN - \2Eo cos ipl]Sip, 
dt 

where N^ and sin(^^ are determined by Eqs. ([23]) . Tak- 
ing into account that As'mipl, < for the both signs of 
A, we find the solution of Eqs. ([5^ in the form 



(32) 



7^. 



SN{t) = SN{0)e~''^(^cosujt- 

-25N{0)e-'<'\A\ . 
o(p(t) = • smujt 



■ sinujt 



1/2 
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This result demonstrates the stabihty of the stationary 
state against fluctuations of the phonon density. The 
character of relaxation to the stationary state turns out 
sign-independent of 6N{0). 

The quantitave character of the solution with the 
damping fluctuations remains for the case S(p{0) ^ 0. 



V. SOLUTION FOR THE CASE OF SEVERAL 
LEVELS 

1 . Let us consider the case when several discrete levels 
lie within the parametric interval. First of all, we analyze 
in detail the evolution to the stationary state for two 
levels at 7 ^ 0. In the case of two levels, taking into 
account double degenaracy of each level in Eq. (|2ip and 
introducing notations Nk^ = ^-fc; — Ni] ipk^ — f^ki = 
ifii] i — 1,2^ we obtain 



dNi 
"dT 

d(pi 
dt 



dN2 

dt 
d(p2 

dt 



-2jNi + 2EoNi cos((5i- 4AA''iiV2 sin((^i- (^2); 

= -2^1 - 2Eo sin ipi - 6AN1 - 8AN2- 

- AAN2 cos(<^i - (P2); 
=-2jN2 + 2EoN2 cos ip2+'iANiN2 sm{tpi~ tp2); 

= -2^2 - '^Ea sin ip2 - 6AN2 - 8^iVi- 

- iANi cos(<^i - (P2); 



(33) 



Here the position of the renormalized level is determined 
by the expression (see (|20)) ) 



^1=^1+ 3ANi + AAN2; Ci = t^i 



Wo 



(34) 



and, correspondently, for ^2 with the replacement 1^2. 
The initial conditions are similar to those in the case of 
a single level, namely, A^i,2(0) ~ 1 smLpi{0) « —£^i/Eo. 

The solution of Eqs. ([55)1 demonstrates an interesting 
asymptotic picture. For any set of the parameters in Eqs. 
(|33p . the number of phonons corresponding to one of two 
levels vanishes. Owing to renormalization determined by 
Eq. (|34p . another level proves to be at the edge of the 
parametric interval with the number of phonons equal to 
from Eq. 



The disappearance of phonons at one of the levels is 
due to the following. This level, as a result of renormal- 
ization, proves to be either beyond the parametric region 
in which pumping is absent, or within the interval be- 
side its edge where the damping exceeds the resonance 
pumping. In the both cases the growth of the phonon 
occupation numbers is suppressed by the damping. 

In most cases we have the characteristic picture repre- 
sented in Fig. 2 when the level 1 closer to the left edge of 
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FIG. 2: Phonon numbers A'^1.2 versus s — 2jt for £,i/Eo 
-0.87, = 0.3 and Eo/j = 5. 

160h 




FIG. 3: Phonon numbers A^i,2 versus s = 2yt for (Eq/^) = 5 
and {£,i/Eo) = -0.87, (^/^o) = -0.3. 



the interval at the initial time moment ("the left level") 
reaches the stationary state {Eo,A > 0). 

It is interesting that, owing to renormalization, the 
level 2 closer to the right edge of the parametric interval 
("the right level") asymptotically proves to be far from 
the edge, = 2.75. 

Figure 3 demonstrates another characteristic situation. 
Now the right level ^2 comes to the stationary state. The 
renormalization, due to N2, shifts the level ^1 to the right 
edge of the parametric interval, £,i/Eq = 0.84, where the 
subsequent evolution of the moderate values of Ni is sup- 
pressed by 7. 




FIG. 4: In the shaded region the right level comes to the 
stationary state. In the rest part of the plane [^1; ^2] the left 
level arrives at this state(i5o/7 ~ 5) 

Figure 4 displays two regions in the plane [^1, ^2] for 
the fixed value (-B0/7) = 5. For the positions ^1, 1^2 
from the shaded region, the right level ^2 arrives at the 
stationary state. For the other positions, the left level ^1 



reaches this state. 

The essential remarks should be added. In princi- 
ple, owing to the renormalization, an outer level with 
|^i|/i?o > 1 (^1 < 0) can be indrawn into the parametric 
interval. It is impossible for any ratio -Bo/7, if ^2 < 0. 
This is already clear from Fig. 4 (with decreasing £'0/7 
the shaded region only enlarges). 




FIG. 5: The shaded region covers the level positions at the 
fixed ratio = 0.4, including the levels with ^i/Eo < 

—1, which are damped in the course of the evolution. The 
stationary state is approached by the state with ^2- 

However, it is possible, when ^2 > 0. Moreover, the 
indrawn level can assymptotically approach the edge of 
the parametrical interval forming the stationary state. 
For any ^2 > there is the limited interval 5c (^2) = 
(|Cic|— £-o)/ii^o, when this event is not realized. At S < 5c 
the indrawn level ousts the level ^2 and reaches the sta- 
tionary position as long as the ratio Eq/j exceeds some 
critical value. The interval Sc increases with decreasing 
^2- The direct numerical calculation gives {6c)max ~ 0.6. 

As an illustration, in Fig. 5 wc plot the plane 
[^i/i?o; £-0/7] and distinguish two regions for the fixed 
value £,2/Eo — 0.4. In the shaded region the positions of 
left levels, including the ones with ^i/Eq < —I, which are 
damped in the course of the evolution, are collected. In 
this case the level with ^2 comes to the stationary state. 
The rest part of the plane is occupied by the level posi- 
tions (including again indrawn levels at ^i/Eq < —I) for 
which the left level ends the the evolution in the station- 
ary state with the macroscopic number of phonons (here 
we imply that ^1 < ^2)- 

Actually, all features discussed here are represented in 
Fig.5. 

3. All qualitative features revealed for two levels are 
reproduced for the case of many levels within the para- 
metric interval. First of all, in the course of the evolution 
only one level proves to be in the stationary state. The 
population of other levels is going to zero in all cases. 
The picture of the evolution as itself can demonstrate dif- 
ferent patterns. In this aspect already the case of three 




FIG. 6: Phonon numbers A'^1,2,3 versus s = 2'yt for {^i/Eo) 
-0.5, = -0.1, i^s/Eo) = 0.1 and (S0/7) = 5. 



levels shown as an illustration in Figs. 6-8 is quite in- 
structive. Really, in Fig. 6 the leftmost level 1 is eventu- 
ally found at the stationary state, while two other levels 
are pushed out of the parametric interval, ^2/Eo = 1.83, 
= 2.03. 

30 H 




FIG. 7: Phonon numbers A/'i,2,3 versus s = 2'yt for (£0/7) = 
3, i^i/Eo) = -0.8, {^2/Eo) = -0.4, (6/So) = 0.1. 

In Fig. 7 the intermediate level 2 comes to the station- 
ary state. Now only the rightmost level is ousted from the 
parametric interval {^^/Eq = 1.89). The leftmost level I 
ends its evolution with zero population at ^i/Eq = 0.9 

At last. Figure 8 displays the interesting situation 
when two levels lie beyond the parametric interval, al- 
though they are sufficiently close to the edge. In this 
case the both levels are indrawn into the interval. The 
leftmost level arrives at the stationary state. 

It should be emphasized that in all cases the stationary 
state arising in the course of the evolution has the number 
of phonons and the phase equal to the values determined 




FIG. 8: Phonon numbers A^i,2.3 versus s — 2'yt for (-E0/7) 
5, (Ci/i?o) = -1.4, (6/i5o) = -1.2, (6/i?o) = 0.7. 



by the expressions p3|) . 

To conclude, we note that the character of the evolu- 
tion described above holds and for A < Q. In this case 
the levels move to the left edge of the parametric interval 
in the course of the evolution. 



VI. PERIODICAL STUCTURE AND SOME 
ESTIMATIONS 

1. The self-consistent temporal evolution of interact- 
ing phonons near the parametric resonance results in the 
formation of the stationary state with the macroscopic 
phonon number of about Eq /A in a single quantum state 
with the common phase for all excitations. The station- 
ary state has one of two energies corresponding to the 
right or left edge of the parametric interval. In fact, the 
stationary state is formed by the phonon pairs which are 
a coherent superposition of two states with the wave vec- 
tors kg and —kg where kg = iOg/co and cOs — li;q/2 ± Eq. 

This state generates the stationary space modulation 
of the atomic condensate. This is easy to reveal deter- 
mining the atomic density as Sn =< x ~^{Pi z)x'{pj z) >, 
where x' has the form (jlOp . Rewriting x' in terms of 
phonon operators and integrating over (Pp, we arrive at 
the expression for the modulation of the ID atomic den- 
sity 



dniz) 2N^ 



{u 



2u N'^ 

l^+vljcos2ksZfiij--^-cos2ksZ, (35) 



where n^^^ — N/L. We suppose that |(5n|/n(-^) ^ 1. The 
phase (/7fc does not enter into this expression owing to the 
accepted relation ifk = f-k- 

It should be emphasized that the static periodical 
structure (|35p is generated by the pair correlations be- 
tween phonons and has the amplitude proportional to 
the number of quasiparticles in the stationary phonon 
state with the macroscopic population. 
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2. Let us make some quantitative estimates. For defi- 
niteness, we consider the quasi ID situation. As is known, 
in this case the chemical potential is equal to /-t = ujj_s 
where s = a^^^n^^\ a^^^ is the 3D scattering length (see 
fl^). The quasi ID situation takes place under condi- 
tion s < 1. Determining the stationary magnitude of the 
number of phonons, we use expressions (13), (16) and 
(|23l) supposing Eq » 7,Cfc- As a result, we find 



N{ jTjs; 

.UJQ 



1] = 



\Slu^\ 



Using this value and definition of the chemical potential, 
we obtain the following expression for the amplitude of 
the periodic term in the longitudinal condensate density 



M 



UJQ 



To estimate the number of the phonon levels iVo within 
the parametric interval 2Eq, we take into account that 
the distance between the nearest levels is equal to 5u)\\ — 

21100/ where co = •\//i/2m is the averaged value of the 
longitudinal sound velocity. Then 



5ljj 



V2jl_ 

TT ^/s 



L 



where — y^l/'muj±. Let us consider a set of parame- 
ters acceptable for experiment 

iV^^lO"^; i = 2-10"^cm; ujjl = 10^ sec~^; 
^ = 0.2: rj=m-\ 

For rubidium atomic gas we arrive at the following esti- 
mates 

s = 0.25; w 10^; w 10"^ No - 1. 

We see that all conditions assumed before are satisfied. 

Thus we can conclude that the state with the macro- 
scopic number of phonons can be realized for a realistic 
set of parameters. 



VII. CONCLUSION 

The analysis presented in the paper has revealed the 
nontrivial picture of the temporal evolution of phonons 
created in a Bose-condensed atomic g result of the 

parametric resonance. The phonon-phonon interaction is 
shown to play a decisive role here. The interaction leads 
to the renormalization of phonon levels which increases 
with growing the number of phonons. It turns out that 
in all cases the evolution stops, when one of phonon lev- 
els within the parametric interval approaches the edge 
of this interval (phonon losses are compensated by the 
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permanent pumping). The macroscopic population of a 
single level is characteristic for this event. At the same 
time the population of all other levels is asymptotically 
going to zero. Actually, the stationary level is filled by 
pairs of correlated phonons with the opposite momenta 
ks,~ks and the common phase. Thus the self-consistent 
nonlinear evolution under consideration has completed 
with the formation of the condensate composed of two 
parts. One part is still in the initial state. The second 
part is made up of the macroscopic number of phonons 
and gives rise to the modulation in the atomic density. 

It should be noted that, for the finite but sufficiently 
low temperatures T < hujQ^ when the initial number of 
phonons with ujk ujq is limited, the results obtained 
remain valid. 

This work is supported by the Russian Foundation 
for Basic Research (Grant No.07-02-00067a) and by the 
Grant for Russian Science Schools (No.NS-6869.2006.2). 



* |manakova@kurm.polyn.kiae.su| 



[1] Yu.Kagan and L.A.Maksimov, Phys.Rev.A, 64, 053610 
(2001). 

[2] Yu.Kagan and L.A.Maksimov, Phys.Lett.A, 317, 477 
(2003). 

[3] F. Chevy, V.Bretin, P.Rosenbusch, K.W.Madison, and 
J.Dalibard, Phys.Rev.Lett., 88, 250402 (2002). 

[4] M.Kraemmer, C.Tozzo, and F.Dalfovo, Phys.Rev.A, 71, 
023613 (2005). 

[5] K. Staliunas, S.Longhi, and G.J.de Valcarcel, 

Phys.Rev.A, 70, 011601(R) (2004). 
[6] M.Modugno, C.Tozzo, and F.Dalfovo, 

cond-mat/0605183 
[7] Yu.Kagan and L.A.Manakova, |cond-mat/0609159 

Phys.Lett.A, 361, 401 (2007). 
[8] Yu.Kagan, E.L.Surkov and G.V.Shlyapnikov, 

Phys.Rev.A, 54, R1753 (1996). 
[9] E.Zaremba, Phys.Rev.A, 57, 518 (1998). 
[10] S.Stringari, Phys.Rev.A, 58, 2385 (1998). 
[11] E.M.Lifshitz and L.P.Pitaevskii, Statistical Physics 

(Pergamon, Oxford, 1980), Part 2. 
[12] M.Wagner, Unitary Transformation in Solid State 

Physics, North-Holland, Amsterdam, 1986 
[13] M. Abromovitz and I.A.Stegun, Handbook of Mathemat- 
ical Functions, Dover, New York, 1970. 
[14] M.Olshanii, Phys.Rev.Lett., 81, 938 (1998). 



